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■ We give a brief review of modern theoretical understanding of the physics of QCD 
at finite temperature and density. We concentrate in particular on the properties of 
hadron gas at relatively small temperatures and discuss in details the physics of the 
phase transition. We argue that the phase transition in temperature is probably absent, 



but it can appear with a vengeance when finite density effects are taken into account. We 
notice also that the temperatures one can expect to reach at the heavy ion collider at 
RHIC are not yet high enough for the perturbation theory in the QCD coupling constant 
to work well. 

43 ■ 1. Introduction. 



The properties of QCD medium at finite temperature have been the subject of intense 
study during the last 15-20 years. It was realized that the properties of the medium 
undergo a drastic change as the temperature increases. At low temperatures, the system 
presents a gas of colorless hadron states: the eigenstates of the QCD hamiltonian at 
zero temperature. When the temperature is small, this gas is composed mainly of pions 
- other mesons and baryons have higher mass and their admixture in the medium is 
exponentially small ~ exp{— M/T}. At small temperature, also the pion density is small 
— the gas is rarefied and pions practically do not interact with each other. 

However, when the temperature increases, the pion density grows, the interaction be- 
comes strong, and also other strongly interacting hadrons appear in the medium. For 
temperatures of order T ~ 150 Mev and higher, the interaction is so strong that the 
hadron states do not present a convenient basis to describes the properties of the medium 
anymore, and no analytic calculation is possible. 

On the other hand, when the temperature is very high, much higher than the charac- 
teristic hadron scale fihadr ~ 0.5 Gev, theoretical analysis becomes possible again. Only 
in this range, the proper basis are not hadron states but quarks and gluons — the ele- 
mentary fields entering the QCD lagrangian. One can say that, at high temperatures, 
hadrons get "ionized" to their basic compounds. In the 0— approximation, the system 
presents the heat bath of freely propagating colored particles. For sure, quarks and gluons 
interact with each other, but at high temperatures the effective coupling constant is small 
a s (T) <C 1 and the effects due to interaction can be taken into account perturbatively. 
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This interaction has the long-distance Coulomb nature, and the properties of the system 
are in many respects very similar to the properties of the usual non-relativistic plasma 
involving charged particles with weak Coulomb interaction. The only difference is that 
quarks and gluons carry not the electric, but color charge. Hence the name: Quark-Gluon 
Plasma (QGP). 

Thus the properties of the system at low and at high temperatures have nothing in 
common. A natural question arises: What is the nature of the transition from low- 
temperature hadron gas to high- temperature quark-gluon plasma ? Is it a phase transition 
? If yes, what is its order ? I want to emphasize that this question is highly non-trivial. 
A drastic change in the properties of the system in a certain temperature range does not 
guarantee the presence of the phase transition point where free energy of the system or 
its specific heat is discontinuous. Recall that there is no phase transition between the 
ordinary gas and the ordinary plasma. 

There are at least 4 reasons why this question is interesting to study: 

1. It is just an amusing theoretical question. 

2. Theoretic conclusions can be checked in lattice numerical experiments. Scores of 
papers devoted to lattice study of thermal properties of QCD have been published. 

3. Some theoretical results can be confronted with the results of the experiments done 
in laboratory with heavy ion colliders. In particular, one could expect to obtain 
interesting results on the collider RHIC which is now under construction. 

4. During the first second of its evolution, our Universe passed through the stage of 
high-T quark-gluon plasma which later cooled down to hadron gas (and eventually to 
dust and stars, of course). It is essential to understand whether the phase transition 
did occur at that time. A strong first-order phase transition would lead to observable 
effects. We know (or almost know — the discussion of this question has not yet 
completely died away) that there were no such transition. But it is important to 
understand why. 

There is also a related question — what are the properties of relatively cold but very 
dense matter and whether there is a phase transition when the chemical potential corre- 
sponding to the baryon charge rather than the temperature is increased. From experimen- 
tal viewpoint this question is no less interesting and important: in the realistic heavy ion 
experiments we have both finite temperature and finite baryon number density. Though 
theoretical predictions are less certain here, we have acquired recently some qualitative 
understanding of high density physics and I am going to discuss it later in my talk. A 
limited time of the talk and a limited volume of this contribution did not allow us to dwell 
upon many interesting aspects of hot QCD physics in details. Also the list of references 
is very incomplete. More detailed discussion can be found in the review 

2. Temperatures are never high. 

The temperature one can hope to achieve in heavy ion collision experiment is not too 
high. For example, the upper bound for the expected temperature in RHIC experiment 
is about Trhic ~ 500 MeV @ and may be twice as much at LHC. 



3 



Let us ask whether the system can already be described in terms of quarks and gluons 
at such temperatures and can be safely called "quark-gluon plasma". In other words: 
whether the effective QCD coupling constant az s (T) is small enough so that the perturba- 
tion theory works well ? We will find out that the answer to this question is unfortunately 
negative (with some reservations). 

Let us consider the basic thermodynamic characteristic of a finite T hadron system - 
its free energy. In the lowest order, we can neglect interaction and the problem is reduced 
to finding the free energy of noninteracting quark and gluon gas. For a pure gluon system, 
the answer can be written immediately using the classical Stephan and Boltzmann result: 
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where (N 2 — 1) is the number of color degrees of freedom. The quarks are fermions and 
integration over momenta with the Fermi rather than Bose distribution gives a somewhat 
different factor. We have 
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where N c Nf is the number of degrees of freedom. Notice that in our World with N c = 
Nf = 3, the lowest order quark contribution to the free energy is roughly 2 times larger 
than the gluon one. In the diagrammatical language, these results correspond to calcu- 
lation of one loop gluon and quark bubbles at finite T. Perturbative corrections to this 
quantity were a subject of intense study during last 20 years [|]. The first correction 
~ a s T A comes from two-loop diagrams and presents no particular complications. It turns 
out, however, that starting from the 3-loop level, we cannot restrict ourselves with a 
given order of perturbation theory but have to sum over a certain class of so called "ring" 
graphs. [] Breaking of naive perturbation series results in the terms in the expansion 
which are not analytic in the coupling constant a s — after the term ~ a s we have the 
term ~ at^' 2 , the next term is ~ a 2 but also ~ a 2 ln(a s ), and the next term is ~ oPj 2 . 
Anyway, the coefficients of all these terms can be found exactly and the result is 



8tt 2 T 4 
45 

F 4 



a a (/i) 
F + F 2 h F 3 



'a s (n) 



3/2 



+ 



7T 



7T 



a. 



71 



5/2 



0{a] lna a ) 



(3) 



where 



1 The reason for that is the bad infrared behavior of the thermal graphs which I have no time to discuss 
here. 
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The expressions ([|) are written for iV c = 3. The coefficients like 237.2 are not a result of 
numerical integration but are expressed via certain special functions. A nice feature of 
the result is its renorm- invariance. The coefficients _F 4 and F 5 involve a logarithmic 
/i-dependence in such a way that the whole sum does not depend on the renormalization 
scale /i. It is important to note that the term ~ o^J 2 in the expansion ([|) is the absolute 
limit beyond which no perturbative calculation is possible. On the level ~ a 3 ., we are 
running into the so called magnetic mass problem, in other words the infrared divergences 
in the graphs become so bad that no recipe for their resummation can be suggested. 

Let us choose \x = 2irT (this is a natural choice, 2irT being the lowest nonzero gluon 
Matsubara frequency ) and Nf = 3. In that case, we have 



1 - 0.9a s + 3.3a 3 J 2 + (7.1 + 3.5 lna,)^ _ 2 0.8a^ /2 
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Note a large numerical coefficient at o> b J 2 . It is rather troublesome because the correction 
~ ct'J 2 overshoots all previous terms up to very high temperatures and, at temperatures 
which can be realistically ever reached at accelerators, makes the whole perturbative 
approach problematic. 

Take T ~ 0.5 GeV. Then 2ttT ~ 3 GeV and a s ~ 0.2. (We use a conservative estimate 
for a s following from T physics. Recent measurements at LEP favor even larger values.). 
The series (ED takes the form 



F = F [l - .18 + .3 + .06 - .37 



(6) 



which is rather unsatisfactory. For T ~ 1 GeV (the upper bound for the expected tem- 
perature at LHC), the situation is just a little bit better. 

We conclude thereby that the hadron matter to be produced at RHIC at LHC cannot 
be called quark-gluon plasma in the technical meaning of this word (a system with weak 
Coulomb-like interaction). It is a kind of "no man's land", a terrain which is difficult 
from the theoretical viewpoint: we cannot describe the system neither in terms of hadrons 
nor in terms of quarks and gluons. Some results can, however, be obtained. First, 
there are numerical results coming from the lattices. Second, some qualitative theoretical 
observations on the nature and characteristics of the phase transition can be made. 
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3. Lukewarm pion gas. 

But before discussing the real experimental situation, consider the system where the 
temperature is small (less than ~ 100 MeV) and the baryon density is zero. Such system 
is well feasible to the exact theoretical analysis. 

Let us first remind the well known facts on the dynamics of QCD at zero temperature. 
Consider YM theory with SU (3) color group and involving Nf massless Dirac fermions 
in the fundamental representation of the group. The fermion part of the lagrangian is 

L f = «E?/7m% (!) 
/ 

where £> M = 8^ — igA a ^t a is the covariant derivative. The lagrangian ([!]) is invariant under 
chiral transformations of fermion fields: 

Q.l,r -> A L:R q L>R , (2) 

where qi,R = |(1 ±7 5 )g are flavor vectors with Nf components and Al,r are two different 
U(Nf) matrices. Thus the symmetry of the classical lagrangian is Ui{Nf) ® Un(Nf). Not 
all Nother currents corresponding to this symmetry are conserved in the full quantum 
theory. It is well known that the divergence of the singlet axial current = Ylf Qfltil 5 Qf 
is nonzero due to anomaly: <9 M j^ ~ g 2 e' J ' I/a ^G'f lu G^p Thus the symmetry of quantum 
theory is SU L (N f ) <g> SU R (N f ) <g> U v (l). It is the experimental fact that (for N f = 2,3, 
at least) this symmetry is broken spontaneously down to Uy(Nf). The order parameter 
of this breaking is the chiral quark condensate SU(Nf) matrix 



< QLfqRf >o (3) 



By a proper chiral transformation it can be brought into diagonal form Sjj/ = — (E/2)5ff 
In the following, the term "quark condensate" will be applied to the scalar positive quan- 
tity S. 

For Nf = 3, the spontaneous breaking of chiral symmetry leads to appearance of the 
octet of pseudoscalar Goldstone states in the spectrum. Of course, the quarks are not 
exactly massless in real QCD, and the mass term is not invariant with respect to the 
symmetry (U). As a result, in real world we have the octet of light (but not massless) 
pseudo-goldstone pseudoscalar states (tt,K,T]). But the small mass of pseudogoldstones 
and the large splitting between the massive states of opposite parity (p/Ax, etc.) indicate 
beyond reasonable doubts that the exact chiral symmetry would be broken spontaneously 
in the massless case. As the masses of the strange and, especially, of u- and d- quarks are 
small 0] , the mass term in the lagrangian can be treated as perturbation. E.g. the pion 
mass satisfies the relation 

FlMl = {m u + m d )T 1 (4) 

(F n = 93 Mev is the pion decay constant) and turns to zero in the chiral limit m u ^ — > 0. 

The presence of light pseudogoldstones in the spectrum is of paramount importance 
for the physics of low temperature phase. When we heat the system a little bit, light 
pseudoscalar states ( pions in the first place) are excited, while p-meson, nucleon and 
other massive degrees of freedom are still frozen. We can study analytically the properties 
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of the system at low temperatures when the medium presents a rarefied weakly interacting 
gas of pions with low energies. Their properties are described by the effective nonlinear 
chiral lagrangian 



where U = exp{2it a (p a /F 7T } is the SU(Nf) matrix (0 a are the pseudogoldsone fields), M. 
is the quark mass matrix and the dots stand for higher derivative terms and the terms 
of higher order in quark masses. When the characteristic energy and the quark masses 
are small, the effects due to these terms are suppressed and a perturbation theory (the 
chiral perturbation theory || ) can be developed. The particular nonlinear form of the 
lagrangian ((^) which involves 4-pion interaction vertices etc. is dictated by the symmetry 
considerations. 

When we switch on the temperature which is small, only the lightest degrees of freedom, 
the pions, are excited. Their density is small and their interaction [described by the 
effective lagrangian @] is weak. 

At T = 0, we habitually describe the system in terms of particles displaying themselves 
as asymptotic states. In the heat bath, the notion of asymptotic states (and the notion of 
S-matrix) makes no sense: due to omnipresent medium the would-be asymptotic states 
scatter on the thermal excitations and do not survive until the collision point. Their 
role is taken over by collective excitations. This notion can be easily visualized when 
recalling the familiar college physics problem of propagation of electromagnetic waves in 
the medium. The frequency and the wave vector of such a wave do not satisfy the vacuum 
relation u = c|k| anymore. A refraction index may appear. In plasmas (also in QGP), 
the dispersion law is modified more essentially and a mass gap develops. The amplitude 
of the classical wave decreases with time due to dissipative effects. That means that 
the wave frequency acquires a negative imaginary part (a positive imaginary part would 
correspond to instability). 

Technically, collective excitations display themselves as poles in the thermal Green's 
functions < tc(x)tt(0) > t , < N(x)N(0) > T , < p(x)p(0) > T , etc. The position of the 
pole o; po i e (k) of the Fourrier image of the Green's function is called the dispersive law of 
the collective excitation with corresponding quantum numbers. 

A problem to find the dispersive law for pion collective excitations is more simple 
theoretically and was solved first ||. As a result, it was found that the pion mass [the 
real part of the pole u>(k — 0)] acquires a temperature correction 



We see that the pion mass slightly increases with temperature, but, if the quark masses 
and M n (0) would be zero, pion mass would remain zero also at finite temperature. This is 
easy to understand: massless pions are Goldstone particles appearing due to spontaneous 
breaking of chiral symmetry. But if the temperature is small, chiral symmetry is still 
spontaneously broken and the massless goldstones should still be there. 




(5) 



Ml{T) = M 2 (0) 1 + 



T 2 
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2 the Green's functions should be retarded ones, but we will not discuss these technicalities here. 



7 



The dispersive laws have not only the real, but also the imaginary part which is physi- 
cally associated with the damping ((k) of the excitations propagating in the medium and 
can be observed experimentally as the broading of corresponding resonance peaks. For 
the pions, damping rapidly grows with temperature. The average value of the damping 
[viz. ((k) integrated over momenta with Bose distribution] was found to be [[jj 

^ T 24F* y> 

Also other hadron states change their properties when temperature is switched on. 
Nucleons, vector mesons etc. are not easily excited when the temperature is small, but 
it makes sense still to study the problem of propagation of a massive state in pion gas 
and find out how the presence of the thermal heat bath affects its dispersive properties. 
Such a study was first carried out for nucleons in ||. The method used was the virial 
expansion over the pion density. In the lowest order, the shift of the pole depends on the 
amplitude of ttN scattering which is very well known from experiment. A brief summary 
of this study is that the nucleon mass practically does not vary with the temperature 
while its damping grows rapidly. At T ~ 150 MeV, the damping is of order ( ~ 100 MeV 
and we cannot reliably calculate and even talk about it at still higher temperatures: the 
description of the system in terms of individual hadrons and also in terms of collective 
excitations with hadron quantum numbers makes no sense in this case. 

4. Vector mesons. CERES experiment. 

Unfortunately, the theoretical results for pions and for nucleons cannot be directly 
confronted with experimental data in heavy ion collisions. 

We hasten to comment that the pole positions of the collective excitations with pion 
and nucleon quantum numbers are quite physical quantities and can be measured in a 
gedanken experiment. 

Suppose we study the spectrum of invariant masses of two 7 emitted from hot hadronic 
matter. The spectrum has a sharp peak associated with the 7r° — > 27 decay. For free 
pions, the width of the peak is very small r° « 8 eV. But thermal pions have a finite 
width Yl = 2 <C >t- Also the maximum of the resonance is shifted towards larger mass 
values according to Eq.(^). Likewise, for nucleons. Proton is believed to decay eventually, 
and most of the decay modes (like e + 7r°) involve hadrons which are stuck within the 
hadronic medium and do not go out undisturbed, but there is also a mode e + 7. For the 
argument sake, let us assume that e + and 7 do not interact with the heat bath, but only 
with the detector. Measuring the distribution in invariant masses of e + 7 and focusing on 
the "nucleon resonance", one could, in principle, extract information on the position of 
the nucleon pole. 

Of course, these gedanken experiments are quite unrealistic. Not even speaking of the 
"nucleon resonance experiment" which is a pure science fiction, also the studying of 27 
spectrum in heavy ion collisions would provide little information about pion properties 
in heat bath. The matter is that 7r° lifetime is at least six orders of magnitude larger 
than the lifetime of the fireball produced in the collision of heavy nuclei . Pions would 
decay on flight in vacuum and their mass and width would be the same as in any other 
experiment. 
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As far as experiment is concerned, the situation is much better for vector mesons, 
especially for the p-meson. Its vacuum lifetime is pretty small, almost all p-s produced in 
the collision decay inside the fireball, and measuring the spectra of e + e~ or p + p~ pairs 
can provide an information about the properties of p in hot hadron medium. Things are 
not so good with u and mesons. Lifetime of hadronic fireball is estimated to be 25 fm/c 
or less ||. It is of the same order as the uj lifetime and almost 2 times less than (f) lifetime. 
A considerable fraction of u-s and most of <p-s would decay outside the fireball. 

Theoretically, the problem to find the dispersive law of the p - meson collective exci- 
tations is somewhat more complicated than for nucleons (the main difficulty is that the 
amplitude of the np scattering in contrast to T n jy amplitude is not known from experi- 
ment), but something can be said, however. First, there is an exact theorem that, as is the 
case for the nucleon mass, the p - meson mass is not shifted to the order ~ T 2 (JTOj] . The 



shift to the order ~ T 4 can be estimated calculating the amplitude T np in the model of 
Ai exchange. The results are similar for nucleon case: the mass is practically not shifted 
with the temperature (at least, while the temperature is not so high that the very notion 
of the p - meson collective excitation loses sense) while the damping (the width of the 
resonance) grows rapidly. (] 

The spectrum of invariant masses of e + e~ pairs in heavy ion collisions has been studied 
experimentally by CERES collaboration. The excess of events with low invariant masses 
which cannot be explained if taking into account only the decay of individual hadron 
resonances has been found. In Fig. [I] we present the brand new CERES plot |12| including 



the separation of data by transverse momenta of the pair. It was found that the whole 
excess is due to the region of small p^ +e . If imposing the cut p^ +e < 500 MeV/c, the 
signal is strong and unequivocal while, with the cut p^ +e > 500 MeV/c, there is no signal 
at all. It was found also that the excess of low mass e + e~ pairs correlates with the square 
of the hot medium density. 

Obviously, this is a collective medium effect which cannot be explained in the model 
with decay of individual resonances. A possible model explanation for that effect is the 
process of 7r + 7r~ annihilation 



7T 7T 



p*(T,p) -> e + e- (8) 



Just the broadening of the p - peak associated with the finite temperature T and the 
finite chemical potential p can explain qualitatively the excess of pairs with low M e + e -, 
its quadratic correlation with the density and the absence of the effect for large p^ +e 
(pions with large Pt are not so abundant in the medium). There is still a problem to 
describe the CERES data quantitatively. May be, under the real experimental conditions 
where the temperature about 150 MeV is achieved, one cannot talk about individual 
hadrons anymore, and CERES result provides us with the first glimps of the no man's 
land between hadron gas and quark-gluon plasma which will be studied extensively at 
RHIC... 



3 The question is still under discussion. In particular, the analysis of Ref.jllJ displays a rapid decrease of 
the mass of p with the temperature. However, this calculation was based on certain model assumptions 
which, in our opinion, do not lead to correct results in this case. 
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Figure 1. CERES spectrum of e + e~ pairs with different transverse momenta. A back- 
ground from individual hadron decays is shown. Solid line is the net contribution of all 
resonances. 



5. Chiral symmetry restoration. 

A physical symmetry which is broken spontaneously at zero temperature is usually 
restored when the temperature is high enough (there are some rather artificial counterex- 
amples which we will not discuss here). Thereby, the chiral symmetry should also be 
restored at high T. At least in massless QCD where the chiral symmetry is the exact 
symmetry of the lagrangian this restoration should be associated with the phase tran- 
sition. We'll come to this soon, but let ask first whether something can be said about 
it by studying the properties of the system in the low temperature phase where chiral 
symmetry is still there and chiral perturbation theory makes sense. 

The important result obtained from the chiral analysis is that the chiral condensate £, 
the order parameter associated with chiral symmetry breaking, decreases with tempera- 
ture. The three-loop CPT calculations have been done in Ref. \TZ . The final result for 
the condensate in massless theory in soft pion approximation (i.e. when the effects due to 
K and r] are disregarded together with effects coming from other resonances) has a rather 
simple form 



T 2 



rp4 



8F2 384Fi 



T 6 , A 
In h 

288F5 T 



(1) 



Here all effects from the higher-derivative term are described by the constant A. Exper- 
imental data on pseudogoldstone interactions give the value A ~ 500 ± 100 Mev. The 
dependence (P together with the curves where only the 1-loop correction oc T 2 or also 




Figure 2. Temperature dependence of the condensate in three loops. 



the 2-loop correction oc T 4 are taken into account (please, do not put attention to the 
"technical" curve marked a° = 0) are drawn in Fig. 0. 

The expansion in the parameter ~ T 2 /8F% makes sense when this parameter is small, 
i.e. when T < 100 — 150 Mev. Strictly speaking, one cannot extrapolate the dependence 
(HD for larger temperatures, especially having in mind that, at T > 150 Mev, the heat 
bath includes a considerable fraction of other than pion hadron states. But as we know 
anyhow that the phase transition associated with restoration of chiral symmetry should 
occur, the estimate of the phase transition temperature (i.e. the temperature when Et 
hits zero, and, as we'll argue a bit later, the transition should be of the second order at 
Nf = 2) based on such an extrapolation is not altogether stupid. This estimate is 

T c « 190 Mev (2) 

or, roughly, one inverse fermi. A more accurate treatment which takes into account 
nonzero m u ^ and also the presence of other mesons in the heat bath gives practically the 
same estimate for the phase transition temperature (assuming, of course, that the phase 
transition is there which, as we will see soon, is probably not the case): nonzero quark 
masses smoothen the temperature dependence of the condensate and push T c up while 
excitation of other degrees of freedom sharpen the temperature dependence and pulls it 
down. These two effects practically cancel each other. 

However, any calculation based on the soft pion technique fails when the temperature 
is of order T c . At such temperatures, we cannot describe the system in terms of goldstone 
fields only (and, indeed, an estimate shows that the contribution of other states in the free 
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energy is not suppressed at T > 150 GeV). Suppose, however, that in the region T ~ T c 
some other effective lagrangian in Ginzburg-Landau spirit can be written which depends 
on the composite colorless fields = qufqif- A general form of the effective potential 
which is invariant under SUi(Nf) ® SU^Nf) is 

V[§\ ~ g{Ti{<&&} + g 2 (Tr{$&}) 2 + g 3 Tr{$&$&} + g 4 (det $ + det $ f ) + . . . (3) 

(the coefficients may be smooth functions of T). Now look at the determinant term. For 
Nf = 2, it is quadratic in fields while, for Nf = 3, it is cubic in fields and the effective 
potential acquires the structure which is characteristic for the systems with first order 
phase transition. A more refined analysis [|14|] shows that the first order phase transition 



is allowed also for Nf > 4, but not for Nf = 2 where the phase transition is of the second 
order. 

Most of the existing lattice data indicate that, at Nf = 3, the phase transition is of the 
first order, indeed, while, at Nf = 2, it is of the second order. 

Up to now, we have concentrated on studying the QCD with massless quarks. But 
the quarks have nonzero masses: m u pa 4 Mev, pa 7 Mev, and m s pa 150 Mev [|J. 
The question arises whether the nonzero masses affect the conclusion on the existence or 
non-existence and the properties of the phase transition. 

Two different experimental (i.e. lattice) works where this question was studied are 
available now, and the results of these two studies drastically disagree with each other. 
The results of Columbia collaboration []nj display high sensitivity of the phase transition 
dynamics to the values of light quark masses. In Fig. [3], a phase diagram of QCD with 
different values of quark masses m s and m u = as drawn in Ref . |T31 is plotted. 



Let us discuss different regions on this plot. When the quark masses are large, quarks 
effectively decouple and we have pure YM theory with SU(3) gauge group where the phase 
transition appears to be of the first order (that is also obtained with some independent 
theoretical arguments). When all the quark masses are zero, the phase transition is also 
of the first order. When masses are shifted from zero a little bit, we still have a first 
order phase transition because a finite discontinuity in energy and other thermodynamic 
quantities cannot disappear at once when external parameters (the quark masses) are 
smoothly changed. 

But when all the masses are nonzero and neither are too small nor too large, phase 
transition is absent. Notice the bold vertical line on the left. When m u = rrid = and m s 
is not too small, we have effectively the theory with two massless quarks and the phase 
transition is of the second order. The experimental values of quark masses (the dashed 
circle in Fig.H) lie close to this line of second order phase transitions but in the region 



where no phase transition occurs. It is the experimental fact as measured in Ref. [15 



12 




Figure 3. Phase diagram of QCD according to Ref.[|l5|]. Mass values (measured in the 
units of the inverse lattice spacing a -1 ) for which the transition is and is not seen are 
denoted respectively by solid circles and squares. Dashed circle corresponds to the physical 
values of masses. 



This statement conforms nicely with a semi-phenomenological theoretical argument of 



Ref. [|Tg] which displays that even if the first order phase transition occurs in QCD, it 
is rather weak. The argument is based on a generalized Clausius-Clapeyron relation. In 
college physics, it is the relation connecting the discontinuity in free energy at the first- 
order phase transition point with the sensitivity of the critical temperature to pressure. 
The Clausius-Clapeyron relation in QCD reads 

disc < qq > Tc = — ^~ dlsc e ( 4 ) 

where disc e is the latent heat. The derivative can the estimated from theoretical and 
experimental information of how other essential properties of QCD depend on m q and 
from the calculation of T - dependence of condensate at low temperature in the framework 
of chiral perturbation theory (see Fig.[| and the discussion thereof). The dependence on 
quark masses is not too weak: dT c /dm q ~ 0.9 — 1.0. From that, assuming that the 



13 



discontinuity in quark condensate is as large as < qq >o, we arrive at an estimate 

disc e < 0.4 GeV/fm 3 (5) 

This upper limit for disc e should be compared with the characteristic energy density of 
the hadron medium at T ~ T c ~ 190 MeV. The latter cannot be calculated analytically 
and we have to rely on numerical estimates. The lattice study in |17| (for the theory with 
2 massless flavors) gives a rather large value e(T ~ 200 MeV) ~ 3 GeV/fm 3 . One can 
safely conclude that the limit (|5|) for the latent heat is several times smaller than e(T c ). 
In reality, the discontinuity in < qq > at the phase transition point (if it is there) is, of 
course, much smaller than the value of the condensate at zero temperature. It is seen from 
the graph in Fig.[2|. The condensate drops significantly still in the region where K and rj 
degrees of freedom are not yet effectively excited. These degrees of freedom (and also p 
and u degrees of freedom) become relevant only at T > 150 MeV when the tendency of the 
condensate to drop down is already well established. Moreover, as was mentioned earlier 
in the paragraph after Eq.(|2|), taking into account these degrees of freedom sharpens the 
condensate dependence rather than smoothens it. 

Thus latent heat of the first order phase transition in the theory with 3 massless quarks 
must be rather small [significantly smaller than the estimate ([5])] which means that the 
phase transition is likely to disappear under a relatively small perturbation due to nonzero 



m s . The numerical findings of Ref.|15| are thereby rather appealing from theoretical 
viewpoint. 

However, as far as numerical lattice experiment is concerned, the question is far from 
being completely resolved by now. A recent lattice study done with Wilson rather 



than Kogut-Susskind fermions displayed quite a different picture shown in Fig.^. 

What is very surprising is the large value of m* ~ 450 MeV (the value of the strange 
quark mass below which the first order phase transition is still possible) as measured in 



Ref . [fig] . If it is true, the point corresponding to our physical world lies well inside the 



first order phase transition region in which case the first order phase transition should be 



of strong variety which contradicts Ref.|15| and also the theoretical arguments above. 

Again, we do not know what of these two measurements is more correct. Possibly, 
intrinsic lattice artifacts are more dangerous when an algorithm with Wilson fermions is 
used but I cannot have an opinion on this point. To resolve the controversy, it would 
make sense to perform a quantitative comparison of the temperature dependence of the 
quark condensate as measured on lattices with the theoretical curve in Fig.|| Such a 
comparison has not been done so far. 

To summarize, at the current level of understanding, the picture where the hadron gas 
goes over to quark-gluon plasma and other way round without any phase transition looks 
more probable. We have instead a sharp crossover in a narrow temperature range which is 
similar in properties to second-order phase transition (the "phase crossover" if you will). 



6. Finite baryon density. Color superconductivity. 

In real experimental situation, not only the temperature, but also the baryon number 
density is high: the colliding nuclei are not only heated, but also squeezed. What do 
the theorists know about the properties of the hadronic system with finite baryon number 
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Figure 4. Phase diagram of QCD according to Ref.|[18||. First order phase transition is 
seen at solid circles and is not seen at blank circles. The star marks out the physical 
values of quark masses. The bend of the second order phase transition line separating the 
region with first order phase transition and the region with no phase transition reflects 
some recent theoretical findings 



density ? First of all, the problem can be solved exactly while the density is still small and 
the medium presents a nuclear liquid. In particular, one can derive the formula describing 



the suppression the quark condensate in dense cold hadronic matter: pO 

2na nN 



< qq > n = < qq > 



FlMl 



+ o{n) 



(6) 



t n ~ 45MeV is the so called sigma-term — the quark contribution to the nucleon 
mass. The correction turns out to be surprisingly high. For the nuclear density n « 
.15(fm)~ 3 , the drop of the condensate constitutes about ~ 30 % ! One should expect 
thereby a drastic change in the properties of the system, probably the phase transition 
in the region n ~ 3rio- Basing on the analogy between (0) and ([]]) and the previous 
discussion for the finite temperature case, one could guess that the phase transition is 
also of the first order. 

This is not true, however, and the phase transition in this case is of the first order. 
Actually, in the imaginary World without electromagnetic interaction, there would be 
two different first order phase transitions in density. The existence of the first such phase 
transition can be qualitatively explained in a very simple way. It is well known that, 
if the protons would not have an electric charge and would not repel each other, the 
most favorable energetically state for the ensemble of a large number of nucleons would 
be the nuclear matter state. It has a definite density n . Suppose we want to study 
the system with zero temperature and an overall density n <uq. The system would then 
present a set of nuclear matter clusters hovering in the empty space. This is what is called 
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the mixed inhomogeneous phase characteristic for the systems with the first order phase 
transition (just imagine a boiling water). The interval of the densities < n < n where 
the system can exist corresponds to one and the same value for the chemical potential 
« c . If we increase the temperature, the nucleons would evaporate from the clusters and 
also pion excitations would show up. While the temperature is small, a system can still 
exist in mixed phase, but it would be less inhomogeneous and the interval of allowed net 
densities An for the mixed phase would be smaller than at T = 0. At a certain critical 
temperature T c , the allowed density interval shrinks to zero and no inhomogenuity in 
nucleon distribution is possible if the temperature is still higher. 

But nuclear matter itself is also an inhomogeneous system: it consists of nucleons where 
the density of matter is higher than in relatively empty regions of space between them. 
Qualitatively, one should expect a first order phase transition associated with squeezing 
the sytem to the point that individual nucleons unite in homogeneous clusters carrying 
large baryon charge. Again, in some density interval the system would exist in a mixed 
phase: dense homogeneous clusters floating in normal nuclear matter or may be clusters 
of the rarified nuclear matter immersed in the dense homogeneous medium, depending 
on what the net density is. Again, the density interval where the mixed phase can exist 
shrinks with temperature and disappears at some critical point. 

Reliable theoretical calculations for the exact position of this second critical point are 
not possible now, but reasonable estimates show that the critical density is ~ (2.5 — 3.0)no- 
Model calculations show that the quark condensate vanishes in the homogeneous phase. 
They also show that, at least in the theory with two massless flavors, it is probably not a 
critical point, but a Uncritical point, i.e. the line of first order phase transitions does not 
end up there but continues as a line of the second order phase transitions going over to 
the finite T second order phase transition at « = 0. 

An exciting recent theoretical finding is that, though the standard quark condensate 
< qq > vanishes in the high density phase, another condensate, the colored diquark 
condensate < qq > which was absent at sero density, appears. This phenomenon has been 
called the color superconductivity. 

A heuristic physical picture for very high densities is rather clear. Remember what 
happens in a usual superconductor. The electrons on the surface of the Fermi sphere 
attract each other due to phonon exchange and form the Cooper pairs. This is associated 
with appearance of the charged order parameter < ee >, Meissner effect, etc. If hadron 
system is very dense, it makes no sense to speak of individual hadrons anymore, the 
system is better described in terms of quarks 0. The colored quarks interact with each 
other. Consider first their perturbative interaction due to 1-gluon exchange. In contrast 
to electrons which always repel each other, quarks repel is their color wave function is 
symmetric (the diquark is in the sextet representation) and attract if it is antisymmetric 
(the anti-triplet representation). Thus, there is a way to form Cooper pairs and they are 
formed. 

That was realized already some time ago [pi], but it was realized also that the pertur- 
bative quark interaction is rather weak. The characteristic values of the superconducting 
gap were estimated to be just several MeV which could not lead to any observable effects. 



4 The temperature is still zero, so we have only quarks and no antiquarks, neither gluons. 
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Figure 5. Phase diagram of hot and dense QCD. 



An important recent observation 22] was that there is also a nonperturbative mechanism 
of forming the Cooper pairs in the attractive anti-triplet channel due to instantons. It 
brings about a rather large values for the mass gap and for the diquark condensate when 
the density is high (above the phase transition), but not asymptotically high. A possible 
form of the net phase diagram for the theory with two massless flavors in the (T, fi) and 
(T, n) planes is shown in Fig.|5| taken from Ref . [ 23| . The solid lines mark out the second 



order phase transitions and the dashed lines — the first order phase transitions. The line 
of the trivial first order phase transition nucleon gas — > nuclear matter is not displayed 
on these plots. The values of T and /x achieved in the real heavy ion collision experiments 
are expected to be close to the position of the upper tricritical point and that may lead 
to some distinct observable effects El. 



7. Disoriented chiral condensate. 

I want to discuss in more details here one particular intriguing possibility that a di- 
rect experimental evidence for the phase transition associated with restoration of chiral 
symmetry can be obtained at RHIC. 

After a head-on collision of two energetic heavy nuclei, a high temperature hadron 
"soup" is created. At RHIC energies, the temperature of the soup would be well above 
the estimate (0) for the phase transition temperature. The high-temperature state created 
in heavy nuclei collision would exist for a very short time, after which it expands, cools 
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Figure 6. Domains of disoriented chiral condensate in cooling hadron soup. 



down and decays eventually into mesons. 

Let us look in more details at the cooling stage. At high temperature, the fermion 
condensate is zero. Below phase transition, it is formed and breaks spontaneously chiral 
symmetry. This breaking means that the vacuum state is not invariant under the chi- 
ral transformations and a direction in isotopic space is distinguished. What particular 
direction — is a matter of chance. This direction is specified by the condensate matrix 
(0). For simplicity, we have assumed up to now that the condensate matrix is diagonal 
Hff = — — 8ffi. But any unitary matrix can be substituted for 5ff (of course, it can be 
brought back in the form 8ffi by a chiral transformation ). In different regions of space, 
cooling occurs independently and flavor directions of condensate are not correlated. As 
a result, domains with different directions of condensate shown in Fig.|6| are formed (cf. 
cooling down of a ferromagnet below the Curie point). 

In our world, we do not observe any domains, however. The direction of the condensate 
in all spatial points is identical. This is a consequence of the fact that u- and d- quarks 
have nonzero masses which break chiral symmetry explicitly, the vacuum energy involves 
a term 

E vac ~ Tr{A^ f S} + c.c. (7) 

and the condensate matrix of the true vacuum state is diagonal. 

However, the masses of u- and d- quarks are rather small and one can expect that the 
domains with "wrong" direction of the condensate are sufficiently developed during the 
cooling stage before they eventually decay into true vacuum with emission of pions. 

This is a crucial assumption. A theoretic estimate of the characteristic size of domains 
they reach before decaying is very difficult and there is no unique opinion on this question 
in the literature. But if this assumption is true, we can expect to observe a very beautiful 
effect p^JT9[| . From the true vacuum viewpoint, a domain with disoriented £/// is a 
classical object — kind of a "soliton" (quotation marks are put here because it is not 
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stable) presenting a coherent superposition of many pions. The mass of this quasi-soliton 
is much larger than the pion mass. The existence of such multi- pion coherent states was 
discussed long ago in pioneer papers [25] but not in relation with thermal phase transition. 

Eventually, these objects decay into pions. Some of the latter are neutral and some 
are charged. As all isotopic orientations of the condensate in the domains are equally 
probable, the average fractions of tc°, tc + , and n~ are equal: < f w o >=< f n ± >= | as is 
also the case for incoherent production of pions in, say, pp collisions where no thermalized 
high-T hadron soup is created. 

But the distribution P(f) over the fraction of, say, neutral pions is quite different 
for coherent vs. incoherent production. In the incoherent case, P(f) is a very narrow 
Poissonic distribution with the central value < f^o >= 1/3. The events with f n o = or 
with f n o = 1 are highly unprobable: -P(O) ~ -P(l) ~ exp{— CN} where N 3> 1 is the 
total number of pions produced. 

For coherent production, the picture is quite different. is proportional to a SU{2) 
matrix. Factorizing over U(l), one can define a unit vector in isotopic space G S 2 . The 
fraction of 7r° produced would be just / = cos 2 6 where 6 is a polar angle on S 2 . The 
probability to have a particular polar angle 6 normalized in the interval < 9 < 7r/2 [ 
the angles 6 > n/2 do not bring about anything new as f(n — 6) = f{6) } is P{0) = sin9. 
After an elementary transformation, we get a normalized probability in terms of /: 

rWf = $j (8) 

As earlier, < / >= 1/3, but the distribution in / is now wide and the values / = and 
/ = 1 are quite probable. 

Thus a hope exists that, in experiments with heavy ion collisions at RHIC, wild fluc- 
tuations in the fractions of neutral and charged pions would be observed. That would 
be a direct experimental indication that a quasi-phase-transit ion occurs where domains 
of disoriented chiral condensate of noticeable size are developed in a cooling stage. One 
can recall in this respect mysterious Centauro events with anomalously large fraction of 
neutral or of charged particles observed in cosmic ray experiments [PEfl . Cosmic ray ex- 
periment is not a heavy ion collision experiment and there are no good theoretical reasons 
to expect the formation of hot hadron matter there, but who knows, may be that still 
was the first experimental observation of the QCD phase transition ? 
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